Turnstile behaviour of the Cooper-pair pump 
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We have experimentally studied the behaviour of the so-caUed Cooper pair pump 
(CPP) with three Josephson junctions, in the hmit of small Josephson coupling 
Ej < Eq. These experiments show that the CPP can be operated as a traditional 
turnstile device yielding a gate-induced current 2ef in the direction of the bias 
voltage, by applying an RF-signal with frequency f to the two gates in phase, while 
residing at the degeneracy node of the gate plane. Accuracy of the CPP during 
this kind of operation was about 3% and the fundamental Landau-Zener limit was 
observed to lie above 20 MHz. We have also measured the current pumped through 
the array by rotating around the degeneracy node in the gate plane. We show 
that this reproduces the turnstile-kind of behavior. To overcome the contradiction 
between the obtained e-periodic DC-modulation and a pure 2e-behaviour in the 
RF-measurements, we base our observations on a general principle that the system 
always minimises its energy. It suggests that if the excess quasiparticles in the system 
have a freedom to tunnel, they will organize themselves to the configuration yielding 
the highest current. 

I. INTRODUCTION 

During the last two decades lots of studies, both theoretical and experimental, have been 
carried out concerning the parametric pumping of charge, an idea originally introduced 
by Thouless in 1982.— This phenomenon is based on the ability of a propagating potential 
well to carry a charge q through a system. This again makes possible controlled pumping, 
by periodically changing the system parameters at a frequency / to induce propagation of 
charge during every cycle. This yields a DC-current I = qf through the system. 

The parametric pumping of charge can be obtained in many different kinds of devices, 
but most of the attention has been directed towards the systems where the charge q passed 
through during each cycle is quantised at a certain number of electrons, q = —ne, where n 
is an integer and the elementary charge of an electron e=|e| is defined as a positive number 
through the paper. This can be realised, e.g., by semiconductor quantum dots by varying 
the height of the tunnelling barriers^ or by one-dimensional ballistic channels in a so-called 
SAW-pump, where the transport is induced by an acoustoelectric wave (SAW).'^'"' The most 
promising candidate so far is the so-called single electron pump consisting of an array of 
three or more mesoscopic metallic tunnel junctions in the Coulomb blockade regime.^^^ Due 
to Coulomb blockade the number of electrons in the islands of the array is very accurately 
controlled and the pumping can be induced by phase-shifted gate voltages,^^'^ which yield 



a current / = —nef. Here / is the frequency of the RF-signal apphed to the gates and the 
integer number n depends on the amphtude of the operating trajectory in the gate variables. 
These devices are more accurate than those based on semiconductors, even sufficient for 
metrological apphcations.-^ In the past the single electron pump has been proposed to close 
the metrological triangle by providing the standard for electrical current 

The only drawback in the single electron pump is the low operating frequency / < 5 MHz, 
which cannot provide high enough current for a current standard. At higher frequencies the 
accuracy is lost due to coherent higher order charge transfer processes known as (in)elastic 
co-tunneling2*i2*ii and other sources of error discussed in Refs. EiEiQIli The maximum 



operating frequency could possibly be pushed higher by using a Cooper pair pump (CPP) 
consisting of three or more small Josephson junctions in series. The pumping in the CPP 
is achieved similarly by gate voltages and the charge is quantised at discrete numbers of 
Cooper pairs in the islands yielding a current / = —n2ef. Due to the different nature of 
Cooper pair tunnelling the operating frequency is now limited by the Landau-Zener (LZ) 
transitions,— which yields several 100 MHz for the upper limit of the operation frequency, 
depending on the parameters of the device. Due to the coherent nature of Cooper pair 
tunnelling, the CPP is also subject to intensive cotunnelling, which reduces the accuracy 
significantly in short arrays.-^^*^ To increase the accuracy one should use longer arrays or 
suppress the cotunnelling by other means.-^^^ 

Another reason for interest in Cooper pair pumping is quantum computing, where Joseph- 
son junction circuits look promising due to their easy scalability and relatively long de- 
coherence times.' 

2iM^ It has been proposed that two capacitively coupled Josephson 
junction arrays could form a quantum bit (qubit) and CPP could also be used to transport 
the information in a more complicated device.-^ When operated in suitable electromagnetic 
environment the CPP could also provide a method to directly measure the decoherence 
rate.A'^'^8- 

In this article we concentrate on transport properties of the CPP and how the non- 
idealities, e.g., dissipation, quasiparticle tunnelling and strong cotunnelhng affect its perfor- 
mance. 



II. COOPER PAIR PUMP 



The Cooper pair pump (CPP) consists of three or more mesoscopic Josephson junctions 
in series with gate voltages capacitively coupled to each of the islands between the junctions. 
It is characterized by several important energy scales. The first one is the typical Josephson 
coupling energy £'j of Cooper pair tunnelling through the junctions and the second one is the 
charging energy Eq due to the small islands between the junctions. The ratio of these two, 
Ej/Ec, is an important parameter in determining the behaviour of the device. To be able to 
obtain pumping of single charges, we have to restrict ourselves to the limit Ej <^ Eq, where 
the dynamics are mainly determined by discrete tunnelling of charge carriers, i.e.. Cooper 
pairs. In addition, the charging energy has to be larger than the thermal energy k-^T to 
prevent thermal excitations, but smaller than the superconducting gap A in the density of 
quasiparticle states to prevent quasiparticle poisoning . This yields the usual chain of 




FIG. 1: (a) A superconducting array of three Josephson junctions (CPP). Here and Ej i^. are 
the capacitance and the Josephson energy of the A:th junction, respectively, (b) StabiUty diagram 
of a uniform CPP at zero bias ^ = on the plane determined by the normalised gate charges 
Qi = Vg,iCg^i/2e. The stable configuration inside each of the hexagons is shown by the kets \nin2)- 



inequalities k-^T < <ti Eq < A. 

In this article we consider a three junction CPP whose circuit schematics are shown in 
Fig. m^a). Each junction has a capacitance Ck and a Josephson energy Ej^k- Gate voltages 
Vg^k are assumed to be externally operated and coupled to islands with capacitances Cg^k- 
The Josephson phase difference ip across the array is related to the bias voltage V according 
to the relation dip/dt = {—2e)V/h. 

Including only Josephson and charging energies and neglecting quasiparticle tunnelling 
as well as other degrees of freedom, the Hamiltonian of the system can be written bsf^^ 

H = Hc^{q,,q2)+Hj, (1) 

where the charging Hamiltonian Hch{<li,<l2) depends on the normalised gate charges, 
Qi = Vg,iCg,i/'2e, and the number of Cooper pairs on each island, rii, according to 
{nin2\Hch{qi,q2)\n[n2)^ = Er.h(uA,U2)S„.. rr[S„ ., „i, where Ui = rn - qi. This model Hamilto- 
nian is described in detail in Refs. lla and|S2- The function i?ch('Ui,M2) gives the details of 
the charging energy and in the presence of the bias voltage V it attains a form 

- 2pel^, (2) 
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where we have assumed a symmetric array with Ci = €2= C^ = C and Cg^i = Cg^2 = Cg- 
Parameter Ag = 2 + Cg/C, Eq = e^/2C is the unit of charging energy, and p is the number 
of Cooper pairs tunneled through the whole array. The Josephson (tunnelling) Hamiltonian 
is given by 

3 

Hj = - ^ Ej^k cos(v9fe), (3) 

k=l 

where Ej^k and <fk are the Josephson coupling energy and the phase difference across the 
fcth junction, respectively. 

At zero bias, = 0, Eq. (0) yields a honeycomb like stability diagram shown in 
Fig. n]^b) .^^i^^ Inside each hexagon the system is stable and there is one charge state \n1n2), 



i.e., the eigenstate of the charging Hamiltonian, Hch, as a ground state. At the edges of 
hexagons and at the triple nodes, two or three charge states are degenerate, respectively. 
The effect of the Josephson coupling, Hj, is to eliminate these degeneracies by introducing 
a coupling between different charge states. This induces new eigenstates which are superpo- 
sitions of charge states and an energy gap will open between different eigenstates near the 
degeneracy lines of the charging part of the energy . 

Tunnelling events in CPP take place as a coherent tunnelling of Cooper pairs, in which the 
system travels adiabatically from the initial charge state to the final one along the eigenstate 
of the full Hamiltonian of Eq. ((T)). This eigenstate is a superposition of the charge states 
with coefficients varying continuously when changing the gate charges, thus resulting in 
a tunnelling of a Cooper pair when a resonance in {qi,q2) plane is passed.—"^ Here the 
resonance means that the initial and the final charge state have the same energy Eq^. These 
resonances for coherent Cooper pair tunnelling are shown as dotted and solid lines in Fig.|21at 
^ = and V > 0, respectively. Similar resonances for the second order process, cotunnelling 
of a Cooper pair, are shown by dashed lines. By cotunneling we mean the coherent tunnelling 
through two junctions simultaneously which is qualitatively similar to the cotunnelling in 
normal state.i^*i2»ii Still higher order processes are weak and play insignificant role with our 
sample parameters. These resonances of coherent tunnelling and cotunnelling overlap in 
case of \^ = 0, but separate as non-zero bias voltage V is applied. Thus, each degeneracy 
node is split into a triangle determined by the first order tunnelling resonances, as shown 
by solid lines in Fig. ISl^--'^^ Inside the triangles the state of the CPP depends on the path 
along which the system has reached the point, thus opening the possibility for hysteretic 
behavior. 

Pumping of Cooper pairs at V = is achieved by adiabatically varying the gate voltages 
along the path encircling one of the degeneracy nodes. The principle of pumping Cooper 
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pairs is explained, e.g., in Refs. |3l[ IitI and Il8l The obtained adiabatic evolution of the 




FIG. 2: Zoomed view of one of the nodes in the stability diagram of the CPP (see Fig. E^b)). The 
dotted lines correspond to 1^ = 0. Solid and dashed lines show the resonance condition for Cooper 
pair tunnelling and cotunnelling in the presence of bias voltage CV/2e = 0.1. 



eigenstates splits the transferred charge into two parts :^ The pumped charge Qp and the 
charge Qs carried by the constantly flowing supercurrent Is- The latter one of these can be 
calculated as the derivative of the dynamical phase r]m = ~ /o (-^m(T)/^)c?r in the state 
|m), Qs/2e = —drjm/dip, while the pumped charge is related to the (/^-derivative of Berry's 
phase^^ 7m = i §{m\dm) attained along the pumping path, Qp/2e = —djm/dipM 

III. SAMPLE FABRICATION 

The sample used in experiments was fabricated by e-beam lithography with a conventional 
self-aligning shadow angle evaporation technique. The two different ex situ process steps 
were performed to make it possible to first evaporate large scale structures containing all 
the contact pads for bias and gate lines and the large guard planes isolating these various 
entries. These guard planes were permanently grounded by bonding them to the bottom 
of the sample holder. The undesirable capacitive crosscouplings between the lines were 
eliminated this way effectively as will be seen later. The design of the large scale structures 
is shown in Fig.Efa). These large structures were fabricated of gold to ensure good electrical 
contact in ultrasonic bonding of wires on the sample stage, and between the two layers 
fabricated ex situ. The so-called quasiparticle traps in biasing lines near the sample were 
also evaporated at the same time.-^*^ These quasiparticle traps consisting of a 2.5 fim x 
0.6 /^m X 60 nm sheet of gold were placed under the biasing line about 4 fim away from the 
outermost junctions. 

The sample itself is surrounded by these large scale structures and it consists of three 
Josephson junctions of about 100 nm x 100 nm in area. The two ~ 1 fim long islands had an 
interdigital type of design to improve the gate coupling and to suppress the crosscouplings. 
The gate lines were also brought to the sample via a 1.8 fim wide channel between the guard 
planes. The small structures were fabricated out of aluminium with aluminium oxide as a 
barrier in the tunnel junctions. 

IV. CHARACTERISATION OF THE SAMPLE 

A. Experimental setup 

The experimental setup is shown in Fig. El^a). All measurements were done in an S.H.E. 
Corporation DRI-420 dilution refrigerator whose minimum temperature is ~ 10 mK. Lower 
parts of the cryostat are surrounded by lead over the vacuum jacket in the helium bath for 
magnetic shielding. 

The fridge has 14 highly filtered lines for DC-signals (Fig. ID (a)). These lines include 3 
stage low pass filtering at different temperatures. At room temperature we used commercial 
low pass TT-filters (Tusonix 4101, -55 dB at 100 MHz) which were connected directly to the 
top of the cryostat. From room temperature down to 600 mK all DC-signals are fed through 
coaxial cables with Nb as an inner conductor and stainless steel as shielding. Between 600 
mK and 60 mK plate each line has 1.5 m of Thermocoax® cable which also forms the next 




FIG. 3: (a) Scheme of the large scale structures of the sample. Gray areas indicate biasing lines and 
the triangular ones gate lines at the top and the bottom. Black areas show the large guard planes to 
suppress parasitic crosscouplings. Grid-like areas in the contact pads are used in ultrasonic bonding, 
and they improve the mechanical strength of the contact. Also the two bigger alignment marks are 
seen inside the two guard planes, (b) An SEM-image of the small structures of the sample used 
in these experiments. The interdigital type design of the gates lowers the crosscoupling between 
them. The large guard planes as well as quasiparticle traps fall outside the image. 

filtering stage (-200 dB at 20 GHz).— At both ends of Tliermocoax® cables tliere are 1 kO, 
resistors in series to improve filtering at low frequencies (< 1 GHz). Tlie last filtering stage 
is at tlie sample stage at tlie base temperature. Tliese filters were commercial stress gauges 
(KYOWA KFG-2-350-D1-23) squeezed between two ground planes forming continuous RC 
strip line filters. Short sample stage wires were made of Cu and were soldered to a printed 
circuit board (PCB). This again was ultrasonically bonded to the DC electrode of the sample 
with aluminum wires. 

Our refrigerator has 4 lines for RF-signals (Fig. E](a)). At room temperature we used 400 
MHz low pass filters (Mini-circuits SBLP-400) and -6 dB fixed attenuators (Inmet). These 
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FIG. 4: (a) Measurement wiring of the dilution refrigerator, (b) A scliematic picture of tlie 
room temperature RF-signal connections. Tlie dotted line shows the part included in the inte- 
grated domestic circuit. The compensation circuit was not needed in the measurements due to low 
crosscoupling in the sample. 



were directly connected to the top of the cryostat and all other room temperature connections 
were made by using flexible SUHNER Sucoflex 104P cables with SMA connectors. From 
room temperature down to 4 K RF-signals are fed through Cu coaxial cables. At low 
temperatures we use BeCu coaxial cables, except between 600 mK and the sample stage 
we use semirigid stainless steel coaxial cables for better thermal isolation. As a whole the 
RF lines have -33 dB attenuation at low temperatures: -20 dB at 4 K, -10 dB at 600 mK 
and -3 dB at sample stage temperature (Inmet fixed attenuators). At the sample stage we 
use MCX connectors and all other connectors are of SMA type. The sample is directly 
ultrasonically bonded by Al bonding wire to a coplanar transmission line mounted on the 
sample stage. When all these lines are connected to the coldest parts of the refrigerator, the 
base temperature is lifted up to ~ 20 mK, as compared to the ~ 10 mK base temperature 
without these lines. 

All DC-voltage measurements were done by using HMS Electronics model 568 low-noise 



preamlifiers and for current measurements we used DL-Instruments 1211 preamplifier. Botli 
amplifiers were powered by battery sources only. Between preamplifier and data aquisition 
(NI PCI-6036E DAQ card) we use a home made battery powered analog optoisolation to 
avoid ground loops and digital noise in our measurements. In Fig. Efb) we show a schematic 
picture of gating signal connections. For this we used a HP8656B signal generator and di- 
vided the RF-signal by using an INMET 6014-2 power divider. One of these signals was fed 
through a phase shifter while the other went directly into the home made circuit schemati- 
cally consisting of two bias-Ts and a possibility to add negative crosscoupling between the 
two signals to compensate the undesirable capacitive coupling between the gates in the GPP. 
The circuit also contained a high quality RF-circuitboard and optoisolated linking to the 
computer which could be used to program all the gains used in bias-Ts, DC-offsets, com- 
pensation and as a main amplification. The control program also included possibility to 
automate the measurements to some extent.— 

B. Current- Voltage dependence 

In Fig. El we present I- V characteristics of the sample, taken with different combinations 
of gate voltages Vg^i and Vg^2- It shows a sharp rise at the beginning of the quasiparticle 
tunnelling branches at bias voltages V ~ ±6AAi/e ~ ±1.2 mV. Also all four major peaks 
in subgap regime corresponding to different possibilities of Josephson-quasiparticle (JQP) 
-cycles are clearly visible.-^i^^iiMi The gate modulation is most pronounced at the gap edge 
as well as in the region of JQP-peaks. From the asymptotic slope of the I-V curve at 
high voltages we obtain the normal state resistance per junction, i?T — 34 kfl, which yields 
E] = AaiRq^^Rt ~ 19 fieV. Here we have assumed that all junctions are identical and 
that the critical current obeys the Ambegaokar-Baratoff relation . Rq = h/{2eY is the 
resistance quantum for Cooper pairs. The value of the charging energy, Eq ^ 129 /leV, was 
obtained using the depth of the dip in the normal state conductance curve taken at 4.2 K.— 
This yields i?j/£'c~0.15. Since the charging energy is well below the superconducting gap 
Eq < Aai ~ 200 /xeV, the parity effect should not be suppressed.— 

The capacitances of the tunnel junctions, Ct ~ 0.62 fF, can be calculated using the Eq 
measured as described above, while the gate capacitances, Cg^j ~ 70 aF {i = 1,2), were 
obtained from the periodicity of the current modulation at high bias. There exists also 
unavoidable crosscouplings, i.e., capacitive coupling of the gate voltage Vg^i to the second 
island, Cu, and vice verse, C21. These parasitic capacitances can be estimated from the 
slope of any known horizontal or vertical periodic structure in the modulation plane similar 
to the one in Fig. IHl where the current is plotted as a function of the two gate voltages 
Vg^i. We found crosscoupling to be C12 ~ C21 ~ 0.16Cg, which is low enough to make it 
unnecessary to apply any active compensation for that. 

In the large scale figure the current around zero bias is not visible but it is seen in the 
blow up. The supercurrent, i.e., the feature around zero bias, depends strongly on the gate 
voltage. It is peaked around V = 0, getting the maximum at |V^|<^£'c/e (see Fig. El). The 
shape of supercurrent also indicates the presence of the electromagnetic environment Z{uj) 
with low impedance of < Re[Z(Lo')] <^ -Rq-— This agrees well with the impedance of the 
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FIG. 5: I-V characteristics of the measured sample, taken at different combinations of gate voltages 
Vg^i and Vg^2- The inset presents the blow up of the I-V curves near zero bias. 



biasing lines, Re[Z(a;)] ~ \/i/c, which is of the order of the vacuum impedance 377 fl. Here i 

and c are the inductance and the capacitance per unit length, respectively, and the resistance 
of the normal metal sections, ~ 1 fi, is negligible. Resonance peaks symmetrically around 
zero bias a.t ±80 /xV, are also visible in this figure, and they will be discussed later. 



To find the correct working point for pumping we mapped out the gate dependence of 
the supercurrent Is{Vg,i,Vg,2) by applying a small bias voltage to maximize Is and mea- 
suring the current while passing systematically all combinations of gate voltages. Since 
Is = —{2e/h){dH/dip) it is clear that it should follow the honeycomb like stability diagram 
of the CPP (See. Fig. E^b)), getting an increase at every degeneracy line and maximum 
value at the triply degenerate nodes. The modulation curve obtained experimentally is pre- 
sented as a contour plot in Fig. IHfa). It does not show the expected pure honeycomb like 
structure, which indicates the existence of non-equilibrium quasiparticles in our system. Yet, 
the pattern cannot be ascribed to pure quasiparticle tunnelling either. We have observed 
this kind of a modulation pattern in most of the (more than ten) measured samples. 

The effect of quasiparticles tunnelling into and out from the islands is to change the gate 
charge by one electron and thus to shift the 2e-periodic stability diagram half a period (i.e., 
e-period) in the corresponding direction. If the rate of these tunnellings is higher than inverse 
of the measuring time, > 100 /is, the current obtained is an average over many 2e-periodic 
honeycombs separated by half a period in the direction of one of the two gate charges. Since 
the shift due to quasiparticle tunnelling is exactly half the period of the stability diagram, 
two shifts in one direction restores the original honeycomb pattern. This feature can be seen 
in the charging Hamiltonian, -ffch; i-e., tunnelling of two quasiparticles through the same 
junction has the same effect on the charging energy as a tunnelling of a Cooper pair. Thus, 
any quasiparticle configuration (nqp^i, riqp 2) can be reduced to {nqp^j = riqp j mod 2}j=i 2 and 
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FIG. 6: Supercurrent Is as a function of two gate voltages Vg^i and Vg^2- The lighter the color 
is, the higher is the current. The structure is composed of four different 2e-periodic honeycombs 
each corresponding to the expected stability diagram of the CPP. These honeycombs are separated 
by half a period due to different number of quasiparticles, which act like discrete variations in 
the gate charge. Lines drawn in the figure present each one of these different honeycombs in the 
case of a symmetric array and they fit the data exactly, with the gate capacitances Cg^i and the 
crosscouplings, Cu and C21, as fitting parameters. A small bias voltage of F ~ 15 fiV was applied 
during the measurement to sit approximately at the maximum of the supercurrent. 

it is enough to consider only four different honeycombs corresponding to, e.g., quasiparticle 
distributions (0,0), (1,0), (0,1) and (1,1). 

As expected, the measured pattern of Fig. IHl composes of four shifted 2e-periodic hon- 
eycombs each corresponding to the stability diagram of the CPP in the absence of quasi- 
particles. These four honeycombs each displaced by half of the 2e-period are illustrated as 
lines. The lines are fitted to the data with the gate capacitances and crosscouplings as fitting 
parameters and they correspond to the symmetric array. In addition to the regular degen- 
eracy nodes the pattern shows some extra peaks which can be explained as intersections of 
different honeycombs. In these points the supercurrent is high in three out of four different 
quasiparticle configurations, yielding a resonance peak. 

In the preceding paragraphs the measured modulation pattern consisting of four shifted 
honeycombs was explained by quasiparticle tunnelling, which was assumed to happen fully 
stochastically in both time and direction. This is the usual assumption in mesoscopic super- 
conducting devices consisting of Josephson junctions.— However, we can also explain these 
four honeycombs by means of quasiparticle tunnelling and by assuming that the system with 



degrees of freedom evolves via states of the lowest energy. In case of a biased array, the 
system can always lower its energy by increasing the number p of Cooper pairs tunneled 
through it. Thus, it will try to maximize the current. If we consider the system as a whole 
including quasiparticles and assume them to have a freedom to tunnel, they will organize 
themselves to the configuration yielding the highest current. Thus, the quasiparticle tun- 
nelling events are still happening stochastically in time but not in direction. However, this 
means that the quasiparticles do not carry the current themselves, which can be the case 
due to Coulomb blockade and the energy gap A. This possibility of quasiparticles to tunnel 
can be clearly observed in Fig. IHl and it can be justified by the same means as in the fully 
stochastic model.-^ 

According to energy-minimisation the system always changes the quasiparticle configu- 
ration to the one corresponding to maximum supercurrent, while varying the voltages Vg^i 
in the measurement of the gate modulation. Thus, this model yields the same combination 
of four honeycombs for the measured current, as we would obtain if the quasiparticle config- 
urations changed fully stochastically. However, in our model the current is higher because 
it is not the time average over many configurations but the largest possible. Unfortunately, 
by measuring the gate modulation it is impossible to resolve between the fully stochastic 
and the energy-minimisation models and find out which one is the valid explanation. Later 
we argue for the latter explanation based on the results of the RF-measurements. 

However, neither of the previous models for quasiparticle tunnelling explain why some- 
times the 2e-periodicity is seen and sometimes not even the sample parameters should yield 
the clear parity effect. To further examine this, we can make a 'worst case' assumption that 
quasiparticles always have freedom to move (due to, e.g., extra subgap quasiparticle states 
in the samples). Then the quasiparticle tunnelling rate F would be essentially determined 
by the charging energy and the general golden rule expression (in the limit Eq 3> k-BT): 
F ^ Ec/{e'^R^) exp(— Ec/^B^), where = Rt/v"^ and r] ~ 10"^ is the relative density of 
quasiparticle states inside the gap .^^i^^ According to our sample parameters and with tem- 
perature of 30 mK, this yields ~ 10^^ hours for the average time between the quasiparticle 
tunnelling events, which is infinite in the time scale of the measurements. Yet, this time de- 
pends strongly on the temperature and if the electronic temperature of the sample would be 
higher, e.g., 300 mK, which could be the case due to inadequate filtering or thermalisation of 
the measurement lines, the time between the quasiparticle tunnelling events would be ~ 500 
ms resulting in e-periodicity in the measurement. This argument yields the same threshold 
temperature of ~ 250 mK for 2e-periodicity as the earlier measurements in Refs. and 
and could explain the lack of the 2e-periodicity in some cases. 

D. Effect of the bias voltage 

Applying bias voltage V changes the stability diagram so that each degeneracy node 
is split into a triangle as discussed in Chapter |H] These triangles are clearly visible at 
the modulation plane measured at \^ = 64 fiV and shown in Fig. [7|^b). Also theoretically 
calculated degeneracy lines giving a resonance condition for Cooper pair tunnelling, are 
drawn in the figure. These lines are calculated using the parameters obtained earlier for 



the measured sample. Figure [7|^a) shows the same modulation plane but at V ^ 0. The 
triangles are formed around the nodes of each honeycomb in Figs.[7|^a) and El and the 'extra' 
nodes formed of intersections of different honeycombs disappear as they should. 
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FIG. 7: Supercurrent Is as a function of two gate voltages Vg^i and Vg^2 with two different bias 
voltages applied: (a) F ~ (Adjusted to the maximum of the supercurrent) and (b) V = 64 ^V. 
The lighter the color is, the higher is the current. In (b) the theoretically calculated resonance 
conditions for coherent tunnelling of Cooper pairs is presented as lines. These lines coincide with 
the pattern of enhanced current in the experimental data. 

With a more detailed inspection one finds the measured current to be slightly increased 
also at resonance lines for cotunnelling, which is not seen at the (current) scale of Fig. [7|^b). 



Good agreement between theoretical resonance lines and experimental data and also the 
reduction of the current inside the triangles proves the measured current to be mainly carried 
by the supercurrent, Is- Thus, quasiparticles present in the system are not acting as major 
carriers of current. The inelastic tunneling of Cooper pairs, where energy is interchanged 
between the tunnelling Cooper pair and other parts of the system, e.g., electromagnetic 
environment or quasiparticles, is also not important. Both these phenomena should increase 
current inside the whole triangle.-^ 

The first resonance peaks in the I-V curve at ^ 80 /iV (see inset in Fig. Ej) can also 
be explained by the effect of the bias to the stability diagram. At 1^ ~ we have the 
highest current flowing at the triple nodes and at the 'extra' nodes, which are also three 
times degenerate. After increasing V the nodes are split and the maximum current appears 
at the doubly degenerate sites, yielding much lower current than at triple nodes of V ^0. 
The new triply degenerate nodes are formed when the triangles situated close to each other 
overlap upon increasing V and the maximum of the current is reached when vertices of the 
neighboring triangles coincide, as is approximately the case in Fig. [7||^b). The closest nodes 
correspond to different quasiparticle configurations, (0, 0) and (1, 1) and are situated along 
the lines Vg^2 = Vg,i + k/2, where k is an integer number, as seen in Figs. Inland [7|^a). Using 
parameters obtained for this sample we estimate the new triple node to be formed between 
neighbouring zero bias nodes at V = 83 /iV, which is exactly the voltage where the highest 
resonance peaks are located in the experimental I-V curve. 

V. OPERATION AS A TURNSTILE 
A. Principle of operation 

A conventional turnstile for electrons or Cooper pairs consists of two arrays (A^ > 2) of 
junctions connected by a common island in between, whose charging and discharging can 
be controlled by a gate.-^^^ The charging sequence is additionally controlled by the applied 
bias voltage, which determines the direction of the obtained DC-current when an AC- voltage 
is applied to the gate at frequency /. Each gate voltage cycle charges and discharges the 
island unidirectionally, thus generating a current I = —ef. 

The CPP can operate as a turnstile because of the hysteretic behaviour within a finite 
bias triangle, opened around node in the modulation plane. The simplest way to describe 
this behaviour is to consider a path in (51,^2) plane with the constraint qi = q2, exiting 
the triangle at both extremes (see Fig. |HI^a)), and to assume that at every resonance, i.e., 
at each degeneracy of the charging Hamiltonian, Hq^, the system is driven to the state 
with lower energy. This simple reasoning alone is enough to explain the turnstile kind of 
behaviour: within every traversal of the path one Cooper pair is transferred through the 
array in the direction of the bias voltage. This principle of operation involves coherent 
tunnelling, cotunnelling and relaxation. It is pictured in Fig.|Hland explained in more detail 
later. 

To obtain this hysteretic behavior one needs dissipation in the system. This can be 
provided, e.g., by the electromagnetic environment or quasiparticles, which both can absorb 



any amount of energy dissipated in the system. Without dissipation the adiabatic passing 
through the degeneracy hne of charging energy Hqt^ will retain the eigenstate of the system, 
i.e., the gap of is not crossed. This will induce coherent tunneling of one Cooper pair 
in the direction which depends whether the system is in the ground or in the excited state. 
The limit for the adiabaticity is determined by the probability of band crossing, the so-called 
Landau-Zener (LZ) tunneling.-ifi This probability can be written a.^^^ 

where q = \/ qf + is the absolute length in the (gi, ^2) plane. In the absence of dissipation 
the probability for band crossing is equal in both directions, i.e., to excite or relax the 
system, but with dissipation the symmetry breaks. The probability to excite the system 
retains the same amplitude, Plz, but the probability for relaxation PrsI = 1 — Plz + Pul is 
significantly increased and reaches almost unity in the case of strong coupling to a dissipative 
element .^i^ 

How this effect of dissipation makes it possible to use CPP as a turnstile is explained in 
Fig. IHl^b), where the energy diagram of the CPP is plotted along the path with time t as a 
parameter \t 1— > (g^, ^2) : Q'i(^) = Q'2(^)] crossing the triangle as shown in Fig. |Hl^a). To be able 
to properly explain the behaviour we have to include the number of Cooper pairs tunnelled 
through the array to our notation: so, instead of \nxn2j we write \nxn2^'p) . 

If the system is initially at the state |00,0) and we start to increase both gates, i.e., 
moving from left to right along the x-axis in Fig. |Hl^b), we come to the point of resonance for 
tunnelling through two junctions simultaneously, i.e., cotunnelling, to the state |01, 0). This 
resonance condition is indicated by a dashed line in Fig. |Hfa). In the case of cotunnelling 
the total coupling for one Cooper pair to tunnel through two junctions, depends on 
the i?j of the junctions and the energy of the intermediate virtual state. Using the second 
order perturbation theory and the value = 50 /xV used in measurements (see next section) 
we obtain Ef jE^ ^ 0.28 yielding the Landau-Zener frequency /lz(-Ej°) ~ 65 MHz. This 
indicates the adiabatic condition to hold during gate excursion and cotunnelling to happen 
with frequencies smaller than 65 MHz. When continuing along the path, the system is 
retained in state |01,0) due to relaxation at resonances crossed. When coming back along 
the same path the system is driven to the state |00, 1) due to a first order resonance, shown 
by a solid line in Fig. IHl^a), yielding coherent tunnelling and the system is kept there by 
relaxation (see Fig. IHl^b)). Then the cycle starts over again with |00, 1) as an initial state. 
Since the situation is fully asymmetric with respect to bias, the operation carries one Cooper 
pair per cycle through the array in the direction of the applied bias voltage V. 

The most significant source of errors in the cycle described above is the inelastic tun- 
nelling of a Cooper pair, which is also induced by the dissipation and the probability of it 
is proportional to the function P{E)^ which gives the probability density of a tunnelling 
Cooper pair to emit (absorb) energy E {—E) to the electromagnetic environment. In addi- 
tion, the presence of quasiparticles increases this probability by providing another way to 
exchange energy. Examples of possible undesirable inelastic tunnelling events are shown in 
Fig.iniby wavy arrows numbered as 1. As one can figure out from the energy diagram, these 




FIG. 8: (a) The path [t i-^ ('7i,'?2) : Qii^) = Q2it)], used in the turnstile measurements (dotted 
hne). It is centered at the degeneracy node of y = and exits the triangle around this node at 
both extremes, (b) Energy diagram of the CPP plotted along the path in (a). The solid lines 
represent the eigenstates of the system and the dashed lines are energies of the pure charge states 
of -ffchj which are indicated at right by kets. The arrows show the method of transferring one 
Cooper pair per cycle through the CPP in the direction of the bias voltage as in a turnstile. The 
dashed arrows correspond to the relaxation into the lower energy eigenstate at around resonances 
while the solid arrays indicate coherent tunnelling or cotunnelling. The thin vertical dotted lines 
are there to clarify corresponding locations between (a) and (b). The diagram has been calculated 
using parameters Ej/Eq = 0.2 and CV/2e = 0.16. 

are not affecting the final outcome of the gate cycle unless there happens several of them 
during one cycle. This would induce transfer of two or more Cooper pairs during that par- 
ticular cycle as in the example shown by thick solid arrows in Fig. IHl Also similar inelastic 
relaxation events which would instantly destroy the outcome of the excursion, are indicated. 
Fortunately though, these relaxation processes are likely to be largely suppressed due to the 
large number of virtual tunnelling events, i.e., higher order of cotunnelling, needed in them. 
These processes are shown by thin wavy arrows striked out and numbered as 2. Thus, the 




FIG. 9: Evolution of the CPP along the gate path of Fig. E^a). The wavy arrows numbered as 1 
show a couple of examples of possible undesirable inelastic tunnellings of a Cooper pair (or of two 
quasiparticles), which is the most significant source of errors in the cycle of Fig. |Htb). However, 
these are not affecting the final outcome of the gate cycle unless there happens several of them 
during one cycle. The thick solid (straight and wavy) arrows show an example of a cycle with 
several inelastic tunnelling events. The charge transferred in that particular cycle is two Cooper 
pairs, i.e., 4e. Striked out arrows numbered as 2 indicate similar inelastic relaxations which would 
instantly induce errors in pumping. These are fortunately suppressed due to the large number of 
intermediate tunnelling events needed. 

system is fairly rigidly 'locked' to transfer only one Cooper pair per cycle. 

B. Measurements 

To test whether the former line of reasoning holds we first measured the gate dependence 
of the current at a small bias voltage, < V < 80 /iV, and 10 MHz sinusoidal signals added 
on top of the DC voltages Vg^i applied to both gates. The AC-signal was fed to both gates at 
the same phase and the amplitude was 1/6 times the 2e-period, thus corresponding to the 
path drawn in Fig. |Hl^a). The obtained current as a function of DC gate voltages is plotted 
in Fig. ^1 The white dotted lines show the structure measured without the RF signals at 

= 0. Data clearly indicate the enhancement of current inside the areas of the triangles as 
predicted in the previous chapter. 

To obtain more quantitative results we froze the DC gate voltages to the values corre- 
sponding to one of the degeneracy nodes and applied a similar sinusoidal signal as before 
at different frequencies. The full I-V dependence was measured instead of sitting at a fixed 
bias voltage to better find the correct bias voltage, < < 83 fiV, where the applied gate 
path would be optimum for turnstile kind of behavior. Some of these measured I-V curves 
are shown in Fig. ^2 

From the I-V curves we can recognize the resonance point at ^ ^ 80 /xV and the optimum 




Vg,1 (Cg/2e) 

FIG. 10: Current with small bias voltage, 0<V <80 /iV, as a function of two DC voltages applied to 
gates V^^i and Vg^2 with 10 MHz in-phase sinusoidal signals added to the two gates. The amplitude 
of the sine was 1/6 x 2e-period. The measured structure of the modulation without RF signals 
and at y ~ is drawn as white lines. 

operation point is the minimum between that and zero voltage, i.e., around y~±50 fiV. If 
we extract the current at V = 50 /iV and plot it against / we obtain the dependence shown 
in Fig. [121 The dashed line shows the ideal I = —2ef dependence which should be obtained 
if the device operates as a turnstile as explained in the preceding section. 

At low / (< 10 MHz) the current increases as J ~ — 2e/ but starts to lack behind at 
around 25 MHz. To find out whether this behaviour is consistent with LZ crossing, we fitted 
the data using 

/ = -2e/[l-PLz(/)]. (5) 

The fitting parameters used were the starting slope (~ 2e), Landau-Zener frequency /lz, 
and the offset in I. This formula did not fit properly the data, mostly due to the nonzero 
asymptotic slope of our data at high frequencies. To take this extra slope into account we 
modified the fitting function to allow for a finite leak current with linear dependence on /. 
The function which we used to fit the data and which is shown in Fig. is of the form 

/ = -2e/gp [1 - exp (-/lz//)] - 2efQi^ + Jq, (6) 

where the fit parameters are /lz, Qp and Ql, which are the charges transferred and leaked 
during one cycle in units of 2e, respectively, and /q, which is the offset in current. 

The existence of this leak current can be physically justified as the result of undesirable 
inelastic tunnelling of Cooper pairs. As explained earlier, these inelastic tunnelling events 
can happen inside the triangle during the operation and their contribution to the resulting 
current is, due to the fact that at the minimum two tunnelling events per cycle are needed to 
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FIG. 11: The I-V curves measured with DC gate voltages tuned to one of the degeneracy nodes 
and in-phase sinusoidal RF signals added to the two gates. Amplitude used was 1/6 x 2e-period 
and frequencies corresponding to I-V curves from lowest to highest absolute current level are 0.5, 
6, 12, 18, 26, 40 and 60 MHz. 
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FIG. 12: Current in a cycle through the degeneracy node (see Fig. lEfa)) as a function of the 
frequency of the RF signals applied in-phase to the two gates. The bias voltage was V = 50 fiY 
and the AG-amplitude was 1/6 x 2e-period. The solid line is the fit by Eq. @ and values for the 
fit parameters are shown in the first line of Table with the exception of Iq . The inset shows the 
same current with the leak current subtracted and the dashed line shows the ideal 2ef dependence 
predicted by the theory. 

affect the outcome, small but nonzero. If we assume the probability for inelastic tunneling 
Pit to be independent of frequency, we obtain an approximate expression for the leak current, 
-^leak = — [2ePj?p -|- (9(Pj^p)]/, which depends linearly on /. The total current is thus of the 
form / = -2e/[l - Plz(/)] + /leak and Ql ~ P^t in Eq. ©. 

Equation (jH)) fits perfectly the data and the values of the fit parameters obtained using it 



are shown in Table HI The analysis of the data was done both at negative and positive bias 
voltages. As explained earlier, the obtained current changes sign with bias voltage. We also 
repeated the measurement in another node with slightly larger amplitude (1/4 x 2e-period) 
of RF signal and these results are also shown in Table HI 

As seen in Table HI the obtained frequency dependence agrees with the theoretical pre- 
diction of — 2e/ at low frequencies, i.e., Qp ~ 1. However, the Landau-Zener frequency 
obtained, /lz ~ 30 MHz, is lower than that estimated (~ 50 MHz) for cotunnelling. This 
discrepancy can originate from the inhomogeneity of the sample. Since Plz is an exponen- 
tial function of E^, it is the smallest Ej^i in the array that determines the LZ threshold. 
The value /lz ~ 30 MHz yields Ej ^ 17 /leV which is a very reasonable value since the 
corresponding number estimated for a symmetric array is 19 /ieV. For the leaked charge we 
obtained Ql ~0.17, which corresponds to a 17 % probability for an extra Cooper pair to 
leak via inelastic tunnelling events shown in Fig. El during a cycle. 



Ampl. (2e) 


Bias 


Qp (2e) 


/lz (MHz) 


Ql (2e) 


1/6 


V>0 


0.985 ± 0.068 


26.2 ±4.4 


0.127 ±0.046 


1/6 


V<0 


0.956 ± 0.080 


29.1 ±6.5 


0.142 ±0.068 


1/4 


V>0 


0.996 ± 0.049 


35.1 ±4.4 


0.132 ±0.031 


1/4 


V<0 


0.976 ±0.100 


24.7 ±5.0 


0.229 ± 0.030 


1/3 


V>0 


4.0 ±2.3 


3.5 ± 1.8 


0.546 ± 0.029 


1/3 


V<0 


3.8 ± 1.3 


4.4 ± 1.4 


0.481 ± 0.027 



TABLE I: Fit parameters using Eq. (jH)) under various experimental conditions. The current scale 
was shifted to set Iq to zero. 

From the preceding experimental results one can conclude that current in the in-phase 
gate cycles follows the relation A/ = —2ef in data measured with amplitudes crossing only 
one triangle on the gate plane. However, this suggests that quasiparticle tunnelling plays a 
minor role when operating in this regime. This is contradicting the general assumption that 
quasiparticles are tunnelling completely stochastically in both time and direction, which 
was one of the explanations for the four honeycombs on the DC-modulation plane in section 
IIVCI This assumption of all stochastic behaviour should prevent the correct operation of 
our in-phase measurement as long as quasiparticle tunnelling is happening at the rate faster 
than the inverse measuring time < (100 /is)~^, which indeed was assumed to explain the 
result of the DC-modulation measurement. 

If the quasiparticle tunnelling is happening less frequently than the frequency of the 
applied gate signals /, a major reduction in current should be observed due to missed cycles. 
Yet, the whole principle of operation breaks down due to the undefined and constantly 
varying trajectory, if these tunnelling events are more frequent than /. Hence, neither of 
these schemes yields A/ = —2ef. 

However, the relatively accurate behaviour in the in-phase measurements can be explained 
using the energy-minimisation model, which is also consistent with the four honeycombs 
obtained in the DC- modulation measurements, as explained earlier in section HVCl As long 
as the amplitude is small enough that the system stays near one degeneracy node and thus at 



the largest possible current given by the different choices of quasiparticle configurations, no 
quasiparticle tunnelling happens and the system is locked to a certain configuration yielding 
accurate operation. But, as soon as we increase the amplitude too much the quasiparticle 
configuration starts to change and thus prevents the accurate transfer. However, if the 
gate trajectory only shortly goes out of the area of the locked configuration, quasiparticle 
tunnelling can be prevented by the operating frequency / being faster than the time of 
quasiparticle tunnelling, and thus the configuration remains locked. 

Increase in the bias voltage (V>80 /xV) has the same effect as increasing the amplitude, 
since the high current areas (edges of the triangles) corresponding to different quasiparticle 
configurations start to overlap making the choice for the optimum quasiparticle configuration 
unclear. At the bias voltages higher than the optimum operation point ~ 50 /zV the 
triangles of the stability diagram spread too much exceeding the amplitude used and the 
accuracy of the turnstile kind of behaviour is diminished until at V>80 fiV it breaks down. 
The effect of this can be seen in the I-V curves shown in Fig. ^2 which also show that some 
frequency dependence is still retained after the resonance point V>80 fiV, but we have no 
clear model to explain this behaviour due to the existence of many overlapping hysteretic 
areas corresponding to different quasiparticle configurations. 




/(MHz) 

FIG. 13: Current at the degeneracy node as a function of the frequency of the applied in-phase 
RF signal of ~ 1/3 x 2e-period amplitude. The circles correspond to data taken at the optimum 
bias voltage 1/ ~ — 50 /iV and the solid line is the fit by Eq. © with parameters shown on the last 
row of Table Q] 

We also performed measurements using two times larger amplitude (~l/3 x 2e-period) 
of the in-phase RF signal. The results obtained at V = —50 /iV are shown in the last two 
lines of Table HI and plotted in Fig. [T31 with a fit by Eq. (jH)) . In the absence of quasiparticle 
tunnelling and thus with a stable 2e-periodic stability diagram, one would expect a twice 
larger Qp (Qp = 2) than in the data presented, e.g., in Fig. [I2l since the trajectory of gates 
crosses two triangles during one cycle. This trajectory is shown by the thick dashed arrow 
and the corresponding 2e-periodic stability diagram by solid lines in Fig. 

As expected, the measurements did not yield such a slope, but about Qp ~ 4 instead. 
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FIG. 14: Trajectory of gates, shown by the thick dashed arrow, in the in-phase measurement using 
the largest amphtude ~l/3 x 2e-period. Thin dotted hues show two honeycomb patterns atV = 
corresponding to quasiparticle configurations (0,0) and (1, 1). Black circle is the operating point 
tuned by DC gate voltages. Solid and dashed lines show the 2e-periodic stability diagrams at 
V = 50 /iV corresponding to quasiparticle configurations (0,0) and (1,1). 

which can be explained with the energy-minimisation model. As seen in Fig. El we already 
move over to the neighbouring triangle corresponding to a different quasiparticle configu- 
ration. Thus the configuration is changed from (0,0) to (1,1) and back during every cycle. 
These two extra quasiparticles responsible for changing the configuration are naturally al- 
ways driven in the direction of the bias voltage thus yielding extra charge of 2e carried 
during every cycle. If we take into account these two quasiparticles transferred during the 
cycle and the three triangles we cross, one corresponding to quasiparticle configuration (1,1) 
and two to (0,0), as shown in Fig. HH we obtain 8e for the total charge carried per cycle. 
This explains the peculiar result of effectively four Cooper pairs transferred per cycle. 

Measurement yields /lz~4 MHz, which is lower than the one, /lz~14 MHz, estimated 
with the smallest Josephson coupling of Ej ^ 17 /ieV, obtained from the small amplitude 
measurements. This deviation is most likely due to the quasiparticle tunnelling involved 
in the cycle. Quasiparticle transfer has an upper limit of < 5 MHz,ii2ii2iiiii^ which together 
with the Landau-Zener limit of Cooper pair transfer (~ 14 MHz) yields a smoother double 
transition starting from <5 MHz. The leak current is much higher compared to the data 
taken with smaller amplitude and interestingly has the dependence of /leak ~ — e/ which 
might correspond to a leak of one quasiparticle during the cycle, but may as well be accidental 
and can also be explained by 50 % probabihty of leakage of a Cooper pair during a cycle. 

VI. PUMPING MEASUREMENTS 

To complement the turnstile-type of in-phase measurements described in the previous 
chapter, we also measured the same sample with 90° phase-shifted RF signals applied to the 



two gates. This provided a circular path around the chosen degeneracy node. The node was 
again found as described earher. At non-zero bias voltage one should traverse around the 
whole triangle to achieve proper pumping. Thus, nonzero bias voltage should allow pumping 
as long as neighbouring triangles do not overlap. However, these measurements did not yield 
as clear results as the in-phase measurements. 
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< 0.0 
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FIG. 15: I-V curves measured at one of the degeneracy nodes with sinusoidal RF signal added to 
the two gates with 90° phase difference. The amplitude of the RF signal was 1/3 x 2e-period and 
the frequencies were 0.1, 0.5, 1, 3, 6, 8, 10, 12, 14, 18, and 22 MHz from the smallest to the highest 
absolute level of current. 

We measured the pumped current at several different nodes and with different amplitudes 
of the RF-signal: 1/9, 1/6, 1/4 and 1/3 times the 2e-period. In Fig.ll5lsome of the I-V curves 
measured using the largest of these amplitudes are shown. They exhibit clear frequency 
dependence which, in turn, varies with the bias voltage. The most surprising detail in these 
I-V curves is that they reproduce the turnstile kind of behavior, so that the pumped current 
is always in the direction of the bias voltage. Thus out-of-phase RF signals increase the bias 
driven current no matter what direction we wind around the node. 

To analyse the 'pumped' current in more detail we examined the current at fixed voltages 
against frequency. Some examples of these plots are shown in Fig. The general behavior 
of the current as a function of the frequency of the out-of-phase RF signals is similar to that 
observed in the in-phase measurements. The current increases first approximately linearly 
but the slope decreases at around 30 MHz, which is again most likely due to Landau-Zener 
tunnelling. After that the behaviour has no general tendency. Curves corresponding to 
different amplitudes or even to same amplitude but different bias voltages differ a lot from 
each other. The only common feature is that the slope tends to increase again. That is why 
the fitting of the formula similar to Eq. ^ does not give a satisfactory result. 

To find the initial slope and an approximate value of the critical frequency /o ~ /lz for 
LZ-crossover, we used an empirical function 




/(/) = -2e/Qp + 2eQ2 [1 - exp (-///o)] / + /q, 



(7) 
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FIG. 16: Current as a function of the pumping frequency. The radius of the circular trajectory 
of gates is 1/4 x 2e-period in (a) and 1/3 x 2e-period in (b). The data in (a) have been taken at 
V = 365 and the dashed hne corresponds to AI = —2ef. Circles and triangles in (b) have 
been taken at V = 35 jj-Y (< 80/uV) and V = 305 /iV (> 80/iV) yielding the initial slope of ~ 4e 
and 3e, respectively. Solid lines are fits by Eq. ((T)), with parameters given in Table ITll 

which represents an exponential crossover from the initial slope —2eQp to another, —2e{Q-p — 
Q2), with the transit frequency /q. The fitting of this function gives more freedom at higher 
frequencies while it still reproduces the initial slope and /q with good accuracy. Figure IT^ 
also shows the function /(/) fitted to the experimental data. The obtained slopes and critical 
frequencies from the fits using Eq. ((7j) are shown in Table |nj 

The resonance point at V ^ 80 fiV, clearly visible also in Fig. is important in the 
analysis of these data as it was in the case of in-phase -measurements. At the two smallest 
amplitudes, 1/9 and 1/6 times 2e-period, frequency dependence of current was seen only at 
the resonance points. The slope obtained at this voltage is quite close to 2e ~ 3.2 x 10~^^ 
A/Hz (Qp f» 1) with both amplitudes. At the amplitude of 1/4 x 2e-period the plateau in 
current was spread over the whole bias range all the way up to the first JQP peak. At bias 
voltages smaller than the one corresponding to the resonance point the supercurrent masked 
the frequency dependence. Supercurrent is very sensitive to the setting of gate voltages, and 
therefore current in this regime is susceptible to even small changes in the trajectory of gates. 



At the bias voltages higher than the resonance point, systematic and rather stable frequency 
dependence was restored yielding the slope of Qp ~ 1. Data fits at this RF-amplitude and 
bias voltage of V = 365 /iV are shown in Fig. ITCT a). The dashed line shows AI = —2ef. 
At the highest amplitude, 1/3 x 2e-period, the frequency dependence of current is very 
distinctive as seen in the I-V curves in Fig. UHl The obtained frequency dependencies split 
into two categories again. At bias voltages smaller than the resonance one ~ 80 /iV the 
initial slope is approximately Qp ~ 2 while at voltages higher than the resonance one this 
slope is smaller with average at around Qp 1.5. Examples of both of these cases are shown 
in Fig. ITTir b). The obtained critical frequency, /o, is always of the same order ~ 30 MHz 
as in the in-phase measurements, thus further suggesting it to originate from Landau-Zener 
events. 



Ampl. (2e) 


V(/xV) 


Qp(2e) 


/o (MHz) 


1/9 


-83 


0.65 ±0.12 


48 ± 16 


1/9 


83 


0.99 ±0.20 


31.5 ± 7.0 


1/5 


-83 


0.97 ±0.18 


26.4 ±4.6 


1/5 


83 


0.97 ±0.34 


26.2 ±9.0 


1/4 


250 ^ 350 


1.00 ±0.26 


24.0 ±7.0 


1/4 


-350 ^ -70 0.886 ± 0.060 


29.0 ±2.4 


1/4 


365 


1.000 ±0.074 


22.1 ± 1.6 


1/4 


-320 


0.934 ±0.086 


23.1 ± 2.3 


1/4 


310 


1.05 ±0.15 


22.9 ±3.4 


1/3 


75 


1.91 ±0.10 


25.1 ± 1.7 


1/3 


55 


1.986 ±0.062 


28.2 ± 1.2 


1/3 


35 


2.130 ±0.060 


32.4 ± 1.4 


1/3 


-75 


1.529 ±0.094 


24.1 ± 1.9 


1/3 


305 


1.53 ±0.10 


16.6 ± 1.1 



TABLE II: Results of the data fits of the out-of-phase measurements. Parameters as in Eq. ((7|). 

It is noteworthy that out-of-phase gate control does not make a distinction whether the 
phase difference is ±90° or —90°. This is the second feature besides the bias dependence 
indicating that we do not pump charges in the traditional sense. Direction of the obtained 
current is determined only by the bias voltage whereas the pumping frequency and radius 
determine the magnitude. This can be explained in approximately the same way as the 
in-phase measurement, where the system always chose the lowest energy state when pass- 
ing a resonance. When trying to pump against the bias voltage the resonance point for 
cotunnelling to happen is reached before the resonance of single tunnelling as can be seen in 
Fig. El It is likely that in our experiment cotunnelling is so strong that, it always happens 
first. This cycle of cotunnelling yields a transfer of two Cooper pairs during one cycle, which 
is twice as much as in pumping in the direction of bias.-^ However, this was not the case in 
our measurements which yielded |AJ| = 2ef in both directions. This discrepancy could be 
due to ~ 33 % suppression of cotunnelling, which sounds reasonable as we can estimate /lz 



for cotunnelling in case of, e.g., the amplitude of 1/4 x 2e-period, to be of the order of few 
MHz, which would suggest rather high suppression. 

One possible explanation for the discovered behaviour arises from the circular pumping 
trajectory. Due to the structure of the stability diagram, Ej° between the states |10) and 
1 01) along the circular path is suppressed by the factor \/5/2 more as compared to other 
cotunnelling events along the trajectory. This directly leads to the reduction of /lz to 
2/5 X /lz and thus to almost doubling of Plz in the middle junction at / ~ /lz ~ 2 MHz. 
This again could prevent the cotunnelling in the middle junction (but not in the other 
two) and thus yield a ~ 33 % suppression of the total cotunnelling probability, which would 
prevent a cycle transferring one Cooper pair against the bias voltage as in the experiments. 



VII. DISCUSSION 

To further discuss the validity of the energy-minimisation model and the model developed 
for the turnstile-kind of behaviour, we briefly review some of the earlier measurements of 
the CPP or similar devices. Somewhat similar turnstile-kind of behaviour in the case of 



a superconducting SET has been reported earlier in Ref. |33|. The principle of operation 
in that particular experiment was based on a similar asymmetric behaviour of the system 
when passing a resonance point.-^ The main difference besides the number of junctions (two 
in that case), was that in that experiment Cr resistors were embedded near the sample to 
provide the required dissipation. These resistors also suppressed any current near the zero 
voltage and thus the current plateaus were only seen at a finite voltage between Ec/2e and 
2A/e. (Our measurement indicates that no additional resistances are necessary to obtain 
this behaviour.) Relatively accurate result of the experiment is consistent with the energy- 
minimisation model, since the gate trajectory they used in Ref. li^ stayed all the time close 
to a resonance point and thus close to the peak in supercurrent. However, this peak in 
supercurrent was suppressed due to the resistive environment, which might have increased 
the probability of quasiparticle tunnelling, yielding a reduced accuracy in 'pumping': This 
agrees with the experiment as they reported current which was inferior to I = —2ef at the 
plateaus observed. 

Our out-of-phase measurement along the circular path around the degeneracy node also 
yielded current — 2e/ but the direction was defined only by the bias voltage as in the in- 
phase measurements. The proper direction of current, depending on the direction of pumping 
and not on the bias voltage, was observed in the measurements of the Cooper pair pump 
by Geerligs et al^ This difference is most probably due to the smaller ratio Ej/Eq ~ 
0.03 in that experiment, which reduced the probability of cotunnelling significantly. This 
suppression is also consistent with a lower Landau- Zener frequency, /lz, obtained in that 
experiment and the lack of the zero bias supercurrent, which prevented the measurement 
of the DC-modulation. Yet, a small step in a pumped current at \^ = remained as a 
consequence of cotunnelling, but the small height of it indicates much higher suppression of 
cotunnelling than in our experiment. Geerligs et al. also state that the step was significantly 
larger in another sample, even reversing the current, which is what we observe. 



The pumping trajectory used in Ref. |3l| was very wide and according to the energy- 



minimisation model the quasiparticle configuration should have changed several times dur- 
ing each cycle. On the contrary the measurement yielded rather accurate behaviour, which 
would indicate no change of the quasiparticle configuration during the operation. This could 
be explained by the fact that quasiparticles had no freedom to tunnel and a certain config- 
uration remained locked all the tirne, thus resulting in pure 2e-periodic stability diagram, 
which is also reported, e.g., in Ref. |^ The failure of measuring the DC modulation in the 
superconducting state in Ref. EH can be explained by the lack of supercurrent (small Ej), 
since there should not be current flowing inside the opened triangles of the stability diagram 
a.t V ^ 0, which was used to map the degeneracy nodes in the normal state. This can be 
verified by the measurement of Fig . W[h). 

In the measurement of Ref. llOl a similar idea of resistive environment as in Ref. was 
used, how in the case of a CPP. The influence of the resistive environment was to suppress the 
otherwise strong cotunnelling of Cooper pairs. The experimental results showed suppression 
of current near the zero bias as in Ref. |33| and thus the current plateaus were again obtained 
only at higher bias voltages. These plateaus were very weak and yielded a behaviour similar 
to our measurement. The value of the pumped current was approximately —2ef but the 
direction was determined just by the bias voltage. This could be due to inability of the Cr 
resistors to prevent the cotunnelling in spite of the small ratio Ej/Eq ~ 0.013. As stated in 
Ref.Ei the quasiparticles were actively present also in this experiment. 



VIII. CONCLUSIONS 

We have shown that current in both the in-phase and out-of-phase gate cycles follows 
the relation A/ = —2ef in data measured with amplitudes crossing or encircling only one 
triangle on the gate plane. These experiments also demonstrate in practise how Cooper 
pair pump could be used as a turnstile with help of dissipation. A model was developed to 
explain the process and experiments clearly demonstrated this behavior quantitatively. 

The failure to observe pure 2e-periodicity in the DC modulation measurement indi- 
cates an active presence of non-equilibrium quasiparticles in the system, as in many earlier 
experiments.— >^i*^2^ However, the agreement between the model developed for the turnstile- 
kind of behaviour and the measurements suggests that quasiparticle tunnelling has a minor 
effect when operating in this regime, which is contradicting the general assumption that 
quasiparticles are tunnelling completely stochastically. To explain this twofold behavior of 
the system, a model was developed, which is based on a very general tendency of a system 
to always strive for minimum energy. It suggests that if the quasiparticles in the system 
have freedom to move, they will organize themselves to the configuration giving the highest 
current, which continuously lowers the energy of the system. The four honeycombs obtained 
in the DC-modulation measurement are consistent with this model and the relatively accu- 
rate behavior in the RF- measurements can be explained with it, too. The model also agrees 
with the earlier measurements of the Cooper pair pump^^^^i and another similar device.— 

We also measured the current pumped through the array by winding around a degen- 
eracy node along a circular path. This showed the correct magnitude in obtained current 
but the direction was not defined by the direction of pumping but it was rather determined 



by the bias voltage thus reproducing the turnstile kind of behavior. This behaviour could 
be explained in general terms, although not as precisely as in the case of in-phase mea- 
surements, by the similar energy minimisation argument, that the system chooses the state 
corresponding to a lower energy when the resonance is passed adiabatically. This does not 
fully explain the observed magnitude of current |/| = 2ef when pumping in the direction 
against the bias voltage. But, it agrees with a strong cotunnelling in our experiment, which 
is due to the large ratio of Ej/Eq ~ 0.15. 

As a final note we state that it is very unlikely for Cooper pair pump as such to be able 
to provide a current standard or otherwise work with high accuracy. The strong cotun- 
nelling and relaxation, among other uncontrollable processes, tend to degrade the pumping 
cycles. However, there might be ways, such as embedding the sample in a highly resistive 
environment^^ or using a combined flux and charge control , ^'^i^" to overcome these difficul- 
ties. Also the use of the CPP to measure, e.g., decoherence time would need much more 
controlled electromagnetic environment to be successful.— 
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